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Abstract. A gravitational lens (GL)-search program, initiated 
in 1990 at the Nordic Optical Telescope (NOT), has revealed 
several possible GL-candidates among a sample of 168 quasars 
(QSOs), chosen from three lists compiled by C. Hazard, D. 
Reimers and I. Surdej, respectively. Some of these candidates, 
selected for having close companions (within 5 arcseconds), 
were imaged in several filters and their colours compared. Low 

I ^ dispersion spectra of the most promising candidates were also 
^ obtained at the NOT and ESO New Technology Telescope 
(NTT). None of these has proved to be strong candidates of 
gravitational lensing effects. We present this new sample of 
QSOs and combine it with previously published optical QSO 
samples in a statistical analysis to yield constraints on flat cos- 

^-H mologies and galaxy velocity dispersions. Finally, by simulat- 
ing larger samples of quasars and gravitational lenses, we dis- 
cuss how the uncertainties affecting our present results would 
be changed. 
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1. Introduction 

One of the (several) aims of GL-surveys is to perform statis- 
tical analysis of the lensing frequency and of the lensed im- 
age configurations (morphology), to derive various fundamen- 
tal parameters for a universe model. In Section 2, we discuss 
and present our sample of QSOs. We first discuss the correc- 
tion for GL incompleteness in the QSO samples and define 
the principles leading to our selection of GL-candidates. Our 
chances of finding GL-systems depend on how the sample of 
QSOs is chosen. Furthermore, the lensing probability is depen- 
dent upon the detectability limits imposed by the instruments 
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and observational conditions at the site. In Section 3, we in- 
vestigate our abilities to detect GL-systems and characterise 
the NOT angular selection function, defining the area (in the 
magnitude difference-angular separation plane) within which 
a pair of point-like images is considered to be a possible GL 
candidate and used in the analysis. In Section 4, we introduce 
the Singular Isothermal Sphere (SIS) galaxy model, employed 
in our statistical model. Furthermore, we also include the im- 
portant magnification bias, which accounts for the increase in 
the observed frequency of lensing. The possible GL-candidates 
found in the NOT sample are presented in Section 5, where we 
also give a brief discussion of each candidate. Finally in Section 
6 we present the total combined sample and the constraints our 
computations impose on galaxy velocity dispersions and on the 
baryonic-matter density parameter Am (see Carroll et al. 1992). 
We also briefly comment on the influence of dark matter (DM) 
on velocity dispersions in gravitational lensing statistics. The 
decrease in the uncertainties of our results when the sample is 
artificially increased are investigated by applying the statistical 
model on synthetic samples. 

2. The QSO sample 

The observed QSO sample (listed in Table 1 and illustrated in 
Fig. 1) was selected from 3 different QSO lists provided by 
C. Hazard, D. Reimers and I. Surdej, respectively. The latter 
one consists of QSOs chosen from the Veron-Cetty & Veron 
(1991) catalogue, while the other two result from independent 
QSO surveys. Some of the QSOs in the two first samples were 
difficult to identify in the images taken, due to crowded CCD- 
fields. To securely identify the QSO in such cases, we have 
obtained finding charts from the digitised Palomar Atlas with a 
magnitude cutoff at about V ~2l. This allowed us to compare 
the observed field with that of the Palomar Atlas and determine 
the presence (or non-presence) of the target QSO. 

Most of the larger optical GL-surveys optimise the ef- 
fectiveness of finding GLs in their surveys by including in- 
trinsically bright QSOs (My < -27). Such QSOs are often 
called highly luminous quasars (HLQs) (Surdej et al. 1987). 
The strong correlation between lensing probability and abso- 
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Fig. 1. Magnitude-redshift distribution of the 168 QSO NOT sample, 
each object indicated by The small dots represent the QSOs from 

the Veron-Cetty & Veron (1991) and Hewitt & Burbidge (1993) cat- 
alogues. The solid lines indicate the absolute magnitudes —27 & —30 
for an Einstein-de Sitter universe model (Q jw = 1 , f2 a = 0). 



itational lensing, because closely separated GL-systems would 
appear elongated under (sub)-critical seeing conditions 

Depending on the seeing conditions, which dictate the res- 
olution and detection of small separation GLs at ground-based 
observatories, this latter selection effect might play an impor- 
tant role in GL-surveys. We will therefore conduct tests through 
simulations in Section 3 to estimate the effects of limited resolu- 
tion on the detectability of candidates and determine an angular 
selection function which describes the validity of our observa- 
tional conclusions. 

Furthermore, a colour criterion used in some QSO surveys 
may bias the selection of QSOs. A QSO superposed on a (lens- 
ing) galaxy is subject to a reddening effect (from the red colour 
of the galaxy) which may endanger a lensed QSO from being 
recognised as a QSO, thereby affecting the completeness of 
the QSO catalogue. This effect was pointed out by Kochanek 
(1991), where as much as 30% of the gravitationally lensed 
QSOs are believed to be excluded from QSO catalogues due to 
such effects. Surdej et al. (1993) corrected for this effect in their 
computations of the lensing frequencies from GL surveys, thus 
multiplying the optical depth by a correction factor Scat = 0.7. 
This factor is included in the present computations, as most 
of the QSOs investigated here are chosen under similar condi- 
tions. This factor will therefore be included in the optical depth, 
deduced in Section 4. 



lute magnitude is discussed in Surdej et al. (1993). Most of the 
QSOs from our sample (marked by * in Fig. 1) tend to lie be- 
neath the line indicating absolute magnitude My = — 27. The 
irregular pattern of the absolute magnitude lines is due to the 
correction for emission lines appearing in the V-filter for various 
redshifts. The emission line correction, K-correction and com- 
putation of absolute magnitudes are discussed in Veron-Cetty 
& V6ron (1991). 

The selection of a flux limited sample of HLQs to max- 
imise the chances of finding GLs is consistent with the previ- 
ously presented samples by Surdej et al. (1993), Crampton et al. 
(1992), Yee et al. (1993) and Maoz et al. (1993) (hereafter the 
SCYM sample). Among the 645 HLQs from this sample, only 
31 were duplicated in the 168 HLQ NOT sample, leaving 137 
unduplicated quasars. The median apparent magnitude of the 
NOT sample is my ~ 17.5, the median redshift is z ~ 2.0 
and the median absolute magnitude is My ~ —27.7 (using 
Om = l.^A = Oand ffo = 50kms-' Mpc"'). 

Among the lens surveys conducted in the previous years, 
most of the objects have been selected from already existing 
QSO catalogues. It has been suggested that many QSO cat- 
alogues are biased against the inclusion of elongated or dif- 
fuse QSOs, because of the selection routines used in search for 
QSOs. According to Kochanek (1991), some GL candidates 
and real GL-systems are lost in surveys searching for GLs as 
a direct consequence of the inherent selection effects in QSO 
catalogues. Elongated QSOs are interesting candidates to grav- 



3. Angular Selection Function 

Several important issues must be taken into consideration when 
we define the angular selection function (ASF) (Surdej et al. 
1993). First we must estimate the limitations on resolution and 
dynamical range (detectable magnitude differences), affecting 
our ability to resolve close components and large magnitude dif- 
ferences, respectively. The selection function will then simply 
describe our observational window in the angular separation- 
magnitude difference plane. We deduce the ASF for our partic- 
ular survey in the following subsections. 

3.1. Small separations 

Fig. 2 illustrates the distribution of the full width at half max- 
imum (FWHM) values for each of the 229 observations of 
the 168 objects. The seeing is distributed in a range between 
0.5-2.0 arc-second, indicating the highly variable image qual- 
ity typical of ground-based observatories. This implies that each 
potential GL candidate has a different detectability range in the 
angular separation-magnitude difference plane. These differ- 
ences obviously affect the possibility of detecting small sep- 
aration GL-candidates (6 < I'.'O). It is therefore necessary to 
weigh each observed image by its image quality (seeing-FWHM 
value). This is done by estimating the maximum detectable 
magnitude difference between two point-like components as a 

' In fact bearing in mind the highly variable seeing and poor res- 
olution of sky-plates, from which the QSOs are often chosen, the 
incompleteness accusations of QSO catalogues seem sensible. 
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Table 1. The complete observational Ust of the 168 QSOs in the NOT sample. See the notes at the end of this table for a description of the 
columns. 



Object 


RA1950 


Deci950 


F 


FWHM 


z 


mv 


Am 


Sep 


GLc 


Not 


HOOOO-263 


000049.45 


-261959.2 


I 


1.11 


4.11 


17.5 









1 


HOOOO+158 


000325.0 


155304 


R 


1.5 


0.45 


16.4 









1 


H0007+1041 


000728.06 


104125 


R 


0.71 


2.2 


18.0* 









2,3 


H0017+154 


001750.5 


152423 


R 


1.3 


2.01 


18.2 









1 


H0029+073 


002943.35 


072200.0 


R 


0.71 


3.26 


18.4 









1 


H0033+156 


003318.2 


153655 


R 


1.5 


1.16 


18.0 









1 


H0035+0725 


003545.99 


072532.6 


R 


0.58 


3.0 


17.8 









3 


H0037+1206 


003724.43 


120654.2 


R 


0.5 


2.5 


18.5* 









3 


H0055-2659 


005532.5 


-265925 


R 


1.1 


3.67 


17.5 






2 


2 


H0105-2634 


010548.19 


-263420.2 


R 


1.45 


3.5 


17.5 






2 


2 


H0131+154 


013129.7 


152744 


R 


0.84 


1.33 


17.3 









2 


HO 143-0050 


014339 


-005038 


R 


0.74 


3.1 


17.0 









1 


H0146+0146 


014644 


014228 


R 


0.93 


2.9 


17.0 









1 


HS0202+1848 


020241.9 


184811 


R 


1.2 


2.71 


17.5 









1 


HS0211+18 


021143.2 


185840 


R 


1.2 


2.48 


18.0 









1 


H0211+1118 


021159.13 


111804.0 


R 


1.07 


2.0 


16.0 









2,3 


H02 14+0820 


021420.81 


082007.1 


R 


1.16 


2.0 


18.0 











TEX0215+165 


021555.3 


163223.8 


V,I,R 


1.13 


1.9 


18.0 









1 


H0216+0803 


021618.1 


080341 


R 


0.66 


2.993 


18.1 











H02 18+0931 


021813.52 


093153.5 


R 


1.13 


2.0 


16.0 









2,3 


HS0227+0558 


022742.0 


055829 


R 


0.71 


2.06 


17.5* 









1 


Q0229+1309 


022902.63 


130941 


V,R,I 


1.2 


2.07 


17.7 


4.0 


6.7 





1 


Q0308+1902 


030851.8 


190225 


V,I,R 


0.86 


2.84 


18.6 


3.1 


4.4 


1 


1 


H0428-1342 


042820 


-134212 


R 


1.26 


3.2 


17.0 









1 


H0432-1246 


043207.2 


-124608 


R 


1.05 


2.8 


17.0 









1 


H0443-1702 


044340.6 


-170243.4 


R 


0.9 


1.6 


17.0 









1 


H0446-1440 


044619 


-144027 


R 


1.42 


1.6 


17.0 









1 


H0449-1325 


044924 


-132532 


R 


1.0 


3.1 


17.0 









1 


H0449-1645 


044959 


-164509 


R 


1.30 


2.6 


17.0 











H0450-1310 


045054 


-131039 


R 


1.07 


2.3 


16.5 









1 


HS0621+6738 


062138.5 


673837 


R 


0.78 


1.57 


18.0* 









1 


HS0624+6907 


062435.2 


690703 


R 


0.83 


0.37 


13.9 


5.8 


8.0 





1 


HS0626+6745 


062655.3 


674552 


R 


0.86 


0.23 


18.0* 


2.5 


10.0 





1 


Q0636+6801 


063647.6 


680127 


R 


0.94 


3.18 


16.5 









1 


Q0640.3+4489 


064018.2 


445339 


V,R 


0.9 


1.805 


18.2 


-3.0 


6.6 





1 


B20650+37 


065035.2 


370927.1 


R 


0.9 


1.982 


18.0 









1 


HS0701+6405 


070130.3 


640545 


R 


0.92 


1.92 


18.0 









1 


HS07 10+6024 


071023.4 


602441 


R 


0.94 


1.77 


18.0 









1 


HS0727+6342 


072746.2 


634218 


R,l 


0.88 


2.38 


18.5 


3.8 


4.7 


1 


1 


HS0734+5956 


073403.0 


595643 


R 


0.94 


1.75 


16.5 









1 


HS0734+6226 


073413.0 


622658 


R 


0.92 


1.08 


18.0 











H0738+3119 


073800.2 


311902.1 


R 


1.14 


0.63 


16.1 


4.4 


6.5 





1 


HS0740+3222 


074022.1 


322219 


R 


1.1 


1.52 


18.0* 









1 


HS0742+3150 


074230.7 


315016 


V,R 


1.04 


0.46 


16.0 









1 


HS0742+3320 


074246.2 


332055 


R 


0.94 


0.61 


17.7 


2.5 


9.9 





1 


HS0743+6601 


074358.6 


660100 


R 


0.9 


2.19 


17.0 











HS0748+3321 


074841.0 


332104.0 


R 


1.12 


1.932 


18.5 











B20751+29 


075151.0 


294952 


R 


1.04 


2.106 


18.5 











HS075 1+6107 


075158.9 


610747 


V,R,I 


1.3 


2.61 


17.5 


3.75 


5.6 







TEX0759+341 


075933.4 


340746 


R 


0.9 


2.44 


18.5 











HS0804+6218 


080401.7 


621826 


R 


0.9 


1.13 


17.5 











OJ508 


080458.3 


495923 


R 


0.9 


1.433 


17.5 











B20808+28 


080832.2 


285402 


R 


0.9 


1.91 


18.0 
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Object 


RAI950 


Dec 1950 


F 


FWHM 


z 


mv 


Am 


Sep 


GLc 


Note 


H0820+1209 


082037 


120928 


R 


0.85 


1.7 


18.0* 









1 


H0821+1022 


082127 


102245 


R 


1.02 


2.185 


18.5* 









1 


H0822+0849 


082217 


084941 


R 


0.85 


1.8 


17.5* 









2 


H0828+1229 


082830.3 


122955 


R 


0.8 


2.8 


18.0* 









1 


H083 1+1248 


083123.11 


124852.5 


R 


0.94 


2.8 


17.8 









2 


H0834+1056 


083405.4 


105645 


R 


0.74 


1.9 


18.0* 









1 


H0834+0937 


083434.4 


093746 


R 


0.8 


2.24 


18.0* 









1 


H0836+1122 


083649.0 


112242 


R 


0.86 


2.67 


18.8 









1 


CS0199 


083807.6 


355542 


R 


1.16 


1.775 


16.0 









1 


HS0839+2858 


083929.3 


285811 


R 


1.27 


1.34 


18.0 









1 


H0841+1256 


084138.9 


125643 


R 


0.73 


2.2 


18.0 









1 


CSO203 


084230.5 


343154 


R 


0.96 


2.126 


17.0 









1 


H0843+1454 


084317.17 


145423.9 


R 


0.86 


2.28 


18.0 









1 


HS0843+2743 


084333.1 


274344 


R 


1.4 


2.03 


18.5 









1 


H0845+1523 


084549.76 


152310.8 


R 


0.72 


2.4 


18.0 









1 


HS0845+3013 


084559.9 


301344 


R 


0.92 


0.67 


17.0* 









2 


H0846+1540 


084620.58 


154039.5 


R 


0.93 


2.9 


18.0 









2,3 


H0846+1517 


084650.95 


151733.2 


R 


0.86 


2.6 


18.0 









1 


RQ0849+2845 


084905.8 


284516 


R 


0.92 


1.27 


20.2 









1 


RQ0849+2820 


084948.7 


282001 


R 


0.94 


0.2 


18.5 









2 


HS0850+2852 


085022.3 


285233 


V,I,R 


1.2 


2.09 


17.5 


4.5 


3.7 


1 


1 


H0850+1755 


085047.2 


175512.2 


R 


0.72 


3.21 


18.0 









1 


H0853+1953 


085335.45 


195312.0 


R 


1.23 


2.8 


18.0 









1 


H0854+1632 


085459.37 


163236.6 


R 


0.88 


2.54 


18.0 









1 


HS0855+2549 


085503.1 


254928 


R 


0.9 


1.80 


18.0 


4.3 


3.3 


1 


1 


H0903+1534 


090305.39 


153450.3 


R 


0.92 


2.8 


18.0 









1 


H0904+2003 


090429.86 


200316.9 


R 


1.0 


2.54 


18.5 









1 


H0905+1507 


090537.49 


150724.7 


R 


0.92 


3.16 


19.0 


2.7 


3.1 


1 


1 


CS0233 


093632.3 


365348 


R 


1.06 


2.03 


17.0 









2 


HS094544630 


094536.1 


463029 


R 


1.0 


0.99 


17.5* 









2 


HS0945+4646 


094555.0 


464614 


R 


1.2 


1.91 


18.0* 









1 


HS0948+4631 


094814.6 


463154 


R 


1.11 


1.8 


18.0* 









1 


HS 1003+4733 


100325.0 


473345 


R 


1.15 


0.27 


17.5* 









2,3 


HS 1004+45 15 


100431.4 


451531 


R 


1.2 


1.3 


17.0* 









1 


HS 1004+4543 


100449.4 


454306 


R 


1.15 


1.66 


17.5* 









1 


CS038 


100904.6 


295648 


R 


1.0 


2.62 


16.0 









2,3 


CS0261 


101653.6 


355654 


R 


1.04 


2.67 


17.0 









1 


Ql 107+487 


110748.2 


484730 


R 


0.77 


3.00 


16.7 









1 


CSO340 


112332.6 


353636 


R 


0.93 


1.285 


17.0 









1 


B3 11 23+395 


112345.8 


393515.1 


V,I,R 


1.06 


1.47 


16.4 


2.5 


5.4 





1 


Ql 123+434 


112349.4 


432607.4 


R 


1.06 


2.01 


18.4 









1 


CBS42 


112839.8 


292054 


R 


1.11 


1.319 


17.0 









1 


US2694 


113417.2 


300808 


V,I,R 


1.14 


1.858 


18.3 






2 


2 


PGl 148+549 


114842.6 


545413 


R 


0.85 


0.969 


16.1 









2 


PCI 158+4635 


115802.9 


463529 


R 


1.2 


4.733 


22.4 


1.5 


8.9 





1 


Q 12076+399 


120738.6 


395617 


R 


1.3 


2.4 


17.5 


4.0 


7.3 





2,3 


Q1210+1731 


121030.4 


173104 


R 


0.87 


2.537 


17.4 









3 


Q1222.7+22.6 


122247.2 


223704 


V,I,R 


1.15 


2.29 


18.0 


0.8 


4.3 


1 


1 


Q1223+1753 


122335.8 


175325 


R 


0.92 


2.918 


18.1 









1 


Q1230+1627 


123039.4 


162726 


R 


0.98 


2.7 


17.8 









1 


CS0151 


123126.9 


292424 


R 


1.2 


2.011 


16.0 









2 


SPl 


135806.0 


390854 


R 


0.75 


3.28 


17.0 









2 


UM627 


135837.0 


000124 


R 


0.98 


1.87 


16.0 









2 


HS1404+32 


140445 


320331 


R 


1.35 


2.37 


17.5* 









2 


CSO609 


140635.7 


491624 


R 


0.95 


2.13 


17.0 









2,3 


Q1412+0925 


141225.6 


092510 


R,I 


1.05 


1.70 


17.4 


2.5 


6.3 
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Object 


RA1950 


■-^ V' V' X ^kJv^ 


F 


FWHM 


2 


77i V 


Am 


Sep 


GLc 


Note 




141322.7 


372015 




1.54 


2.36 


18.0 









1 


HS1420+S26 


142021 1 

X ^^\J^ X • X 


323647 


R 


0.94 


0.69 


17.5 









2,3 




1421 17.5 


330555 


JJ^ 


1.03 


1 


1 6 5 

i u — ' 








u 


9 3 


1 l_//V 1 *+Z, i. T J J 7 


142144.1 


355606 


]^ 


Q8 


1 57 


17 5 







u 


2 




142947.2 


1 1 5306 


R 


1.85 


3.01 


18.6 













143 SI 3 9 


-01 341 3 




1.7 


1 31 


1 6 
i u.u 











01449+9QS 


144244.4 


293142 


VR I 


1.7 


2.67 


16.2 


9 


3 4 


1 






1 S044S 


542324 


VI R 


84 


1 
i .y 


17.0 


2.4 


6 4 



u 






151708 9 


935659 




06 


1 
i.y 


17 3 







u 




SP43 


1 52030 


412212 


R 


0.84 


3.1 


17.5 











LB9707 


i JZ, JwO .0 


212436 


R 


1.1 


1.024 


18.0 








u 




TFX1 558+1 87 


1 55809 Q 


1 84654 


R 


0.9 


2.4 


18.0 











PKS1 559+17317 


1 ^5904 6 


172236 


R 


1.05 


1.044 


18.0 











JJ J ± \JZ< A. T J 


169193 4 

X \JZ< X Z<^ 


391627 


VI R 


0.6 


1.07 


17.5 


2.3 


3.8 


1 






169345 4 


965354 


R 


1.15 


2.52 


16.0 











HS 1 696+64 


1 69690 4 


643339 

VJT J J JZ. 


R 


0.92 


2.32 


16.5* 


3.2 


9.0 







Ol 698+3808 


1 69898 

X VJZjOZjO.VJ 


380447 


R 


0.73 


1.461 


16.8 











0163fl+374Q 


1 63039 9 

X \J JZf .Zf 


373925 


R 


0.73 


2.037 


18.2 











01633 5+96 7 


1 63334 8 


264417 


R I 


1.16 


1.84 


17 













1fi^4S1 7 


703737 


R 


2.24 


1 33 


14 






n 

u 






i U J Ot-U .0 


612153 


R 


88 


0.46 


1 8 0* 
i o.u 







u 




01643+401 


164396 7 


400443 


R 


0.84 


1.877 


18.0 











HS 1700+64 16 


1 70040 5 


641625 


R 


0.8 


2.74 


16.1 






2 




HS1 704+6048 




604831 


R 


1.1 


0.37 


1 5 3 







u 




1 F 1 704+7 1 


1 70500 6 


7101340 

/ X V/ X .^y ^ . V/ 


R 


1.0 


2.00 


17.5 











HS170Q+6097 


170951.4 


602723 


R 


1.1 


1.54 


17.0* 









. 




171'i30 7 


533194 


VI R Z B 


79 
u. / z. 


1 03 
I.y J 


16 5 


116 

U. 1 iU 


3 7 
J, / 


1 




RS1716+683Q 


171627.9 


683948 


R 


1.7 


0.78 


18.5 











POI 71 8+481 

1 V.J X / X Olio X 


171817.7 


48071 1 


R 


1.45 


1.08 


14.7 











01 790+3470 


1 79009 3 

X / Z.V'V'Z.. J 


344149 


VI R 


1.5 


2.4 


18.4 


1.3 


7.7 







01 799+34 


1 79959 8 

X / Z.Z.-J 7 .0 


344437 


R 


1.2 


2.2 


18.5 


1.9 


6.1 







HS1793+f;'i'iO 


179307 9 


655096 


R 


1 1 5 


1 45 


17 0* 

1 / .u 







u 




TFX1 796+344 


179601 3 

X / Z<\JW X .J 


342504 


I 


1.6 


2.425 


18.5 











HS1739+6535 


1 73946 9 

X / JZ.T^V/.Z. 


fiS'^S?'^ 

vJ J J JZ, J .\J 


R 


0.97 


0.86 


17.6 









. 


HS1749+fS147 


1 7499 1 6 


614711.0 


R 


1.22 


59 

U — 'Z, 


1 8 6 
i o.u 







u 




HS174Q+7nnfS 


1 74Q03 4 


700630 


R 


04 


77 
u. / / 


17.4 







u 




HSl 803+6737 


1 80337 4 


673754 


R 


1.9 


0.14 


15.8 











HS 1807+6048 


180718 5 

X 0\J 1 X O..J 


694857 1 

\jy-rOJ / . X 


R 


1.8 


0.051 


14.4 











PKS1 891 +1 

L XViJ X Oi. X T X vy 


182141.6 


104244 


R 


1.03 


1.36 


17.3 


3.4 


8.0 







HS1 894+6507 


1 89435 9 


650737 


VB R I 


1.75 


0.3 


16.5* 


1.5 


5.2 







4r5fi 98 


1 85730 1 


564146 


VZ R I 


1.2 


1 505 


17 3 


9 7 


3 9 


1 




HS1Q46+7658 


194641.0 


765896 


VI R 


1.2 


3.02 


15.9 











Pr9047+0193 


904750 7 


019356 

VJ X Z< J JVJ 


R 


1.10 


3.700 


19.7 











TFX?n4R+1 Qf^ 

1 Xj^vZ-Ut^OT i 7U 


9n4RS6 7 




R 


85 

U.O J 


2.36 


1 8 5 
i — J 








u 




TFX91 16+903 


91 1617 ? 

Z< X X VJ X / .z< 


909199 

Z<VJZ< X z<z< 


VR I 


0.98 


1.68 


17.0 


2.75 


3.7 


1 




TFX9 197+1 76 


212717.7 


174136 


VR Z I 


0.9 


2.01 


18.0 


1.3 


2.1 


1 




TFX91 97+^48 

1 1-^ /\ _ i Z, / T JT-O 


9 1 9746 9 


345394 


R I 


1 6 
i .u 


2.4 


1 8 5 
i — J 


9 7 


8 1 
0.1 



u 




PKS91 50+050 


91 5054 1 

Z. X -J VJ-Jt. X 


059909 

VJ-JZ.Z.VJ7 


R 


0.9 


1.070 


17.8 











H9930+1 14 

xxzz.JvyTX X 


993007 81 

zz. jvjvj / .0 X 


1 12822.5 


R 


0.75 


1.04 


16.5 











Q2231+0125 


223142.1 


012541 


R 


0.9 


1.0 


18.2 











H2251+158 


225129.4 


155254.0 


R,I 


0.7 


0.86 


15.3 











H2251 + 113 


225140.7 


112036.8 


R 


1.51 


0.32 


15.8 











TEX2300+345 


230059.7 


343038 


R 


0.7 


2.400 


18.0 











B223 10+383 


231036.3 


383123.0 


R 


1.6 


2.17 


17.5 











H23 10+00 18 


231050.8 


001826.4 


R 


0.83 


2.2 


17.0* 
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Object 


RAl95() 


Dec 1950 


F 


I'WHM 


z 


mv 


Am Sep (;Lc 


Note 


H23 14+160 


231409.7 


160143.5 


R 


1.5 


0.66 


18.0 





2 


H2355+0845 


235521.6 


084521.6 


R 


0.7 


2.8 


18.0* 


2 


2 


Q2359+0653 


235906.77 


065313.3 


R 


0.45 


3.24 


18.8 





2,3 


H2359+1305 


235957.4 


130532 


R 


0.77 


2.8 


18.0* 





3 



Notes : 

Col. 1: QSO identification, where the id indicates list of origin; H - from Hazard, HS - from Reimers (Hamburg Stemwarte) and others from 
the Veron catalogue 

Col. 2: QSO RA (hhmmss) in equatorial 1950 coord. 

Col. 3: QSO Dec (ddmmss) in equatorial 1950 coord. 

Col. 4: Filter types 

Col. 5: Seeing (FWHM) in arcsec 

Col. 6: Redshift z of the QSO 

Col. 7: V-band magnitude of the QSO as estimated from our data or tabulated in other QSO catalogues. The QSOs with no previously measured 
magnitude were estimated by comparing them to QSOs on other CCD frames (for a similar exposure time). These estimated magnitudes are 

marked by an asterix (*) and have errors up to ±0.5 mag 

Col. 8: Magnitude difference between closest companion with weaker apparent brightness and within 5" of QSO 
Col. 9: Separation (in arc second) of closest companion with weaker apparent brightness and within 5" of QSO 

Col. 10: A "I" indicates that the QSO is treated as a possible GL candidate, while a "2" indicates that the object is elongated, diffuse or otherwise 
especially interesting. The latter objects are not treated as GL candidates in the present analysis, but should be further analysed in future work 
Col. 11: A "3" indicates that there are no PSFs other than the QSO in the image, a "2" indicates that the QSO appears elongated (this is 
sometimes due to bad tracking) and a "1" indicates that there are one or more PSFs and the QSO is not elongated 



100 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

FWHM 



Fig. 2. The seeing (FWHM in arc-sec.) distribution of 229 observations 
of 168 QSOs 



function of their angular separation, which we evaluate by per- 
forming a simulation-test. Using a typical point spread function 
(PSF) for different seeing values, we have produced simulated 
images of various closely separated pair-scenarios, adopting a 
random flux ratio. The results of these tests are illustrated in 
Fig. 3, which also show the selection function for some of the 
other surveys being added in Section 6. We discuss these tests 
in more detail in Section 3.3. 



3.2. Large separations 

In the identification of possible GL-candidates we chose all 
QSOs with a companion weaker than the QSO within the area 
7rr^. If more than one close companion exists, we chose the 
brighter one, because it has the highest probability of being a 
gravitationally lensed image of the QSO. We have adopted a 
cutoff radius rc =5 arc-second to secure the inclusion of most 
GL-systems produced by a single isolated galaxy and exclude 
those produced by clusters of galaxies. Using a larger cutoff 
radius would also include many star or galaxy companions 
and, hence, increase the number of false positives (Kochanek 
1993b). 

3.3. The ASF 

To determine the ASF of the NOT sample for a wide range of 
seeing values we performed detection tests on synthetic data. 
These data were produced from numerical PSFs obtained from 
one or several stars in an image of a certain quality (we used 
PSFs with FWHM of (y.'6, a'75, 1'.'Oand 1'.'2). The synthetic im- 
ages consisted of a closely separated image-pair of PSFs, which 
were produced using random angular-separation and magnitude 
difference. Photon (Poisson) noise and read-out noise were in- 
cluded in these simulations. Inspection by eye was performed 
for a series of synthetic images in each FWHM category. For 
each of these categories, a set of constants (Ci, C2, A and B) 
were determined, describing the particular ASF for that partic- 
ular seeing (FWHM) value. Cubic interpolation was applied to 
determine these constants for the intermediate seeing values, 
thereby estimating the ASF for any given seeing value within 
reasonable range of the actual measured values. The ASF is 
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valid for seeing values from these tests in the range O'.'S < 
FWHM < I'.'S. We also found from these tests that the detec- 
tion probabilities are slightly dependent on the position angle 
(PA) of the image pair for small separations. Specifically, align- 
ments along the rows or columns of the detector decrease the 
probability of separating the pair on the basis of inspection by 
eye. 

The complete ASF, which we apply in our work, is now 
represented by 



C no solution for 0" < 6» < Ci 
Am(6l) = Ib + AxO for Ci < 61 < C2 
I 5.0 mag for C2 < < Vc 



(1) 



where C\ and C2 are determined for each FWHM value and 
constitute the starting value and "break" point of the ASF, re- 
spectively. Some of the constants that determine the ASF for 
various values of seeing (FWHM) are listed in Table 2. The 
ASF is also shown in Fig. 3. 



Table 2. Constants describing the ASF for various seeing (FWHM) 
values. 



FWHM 


A 


B 


Ci 


C2 


O'.'S 


22.3 IS 


-4.305 


0.192 


0.416 


0'.'6 


17.700 


-3.520 


0.198 


0.482 


0'.'7 


14.193 


-3.185 


0.224 


0.576 


O'.'S 


12.845 


-3.536 


0.276 


0.664 


0'.'9 


13.390 


-4.508 


0.336 


0.710 


I'.'O 


14.176 


-5.712 


0.402 


0.756 


I'.'l 


13.809 


-6.854 


0.496 


0.858 


1'.'2 


12.599 


-8.000 


0.636 


1.032 


1'.'3 


11.526 


-9.3S5 


0.S14 


1.24S 



4. Statistical model 

Following previous GL statistical studies such as Turner 
et al. (1984),hereafterTOG;Fukugita& Turner (1 99 l),hereafter 
F&T;Surdej et al. (1993); Kochanek (1993c) and Maoz & Rix 
( 1 993), we have used the well known singular isothermal sphere 
(SIS) model in the statistical analysis of the combined SCYM 
and NOT samples. In this section we outline the SIS model as 
well as the magnification bias. 

A SIS galaxy model is characterised by its constant deflec- 
tion angle (independent of the impact parameter) giving the 
critical (Einstein) radius bcr = A-Ka^ / c^{Dis / D g), where Dig 
and Dg are the proper distances ^ lens-source and observer- 
source, respectively. A "strong lensing" event occurs when the 
QSO is projected within the area defined by the critical radius, 
7r6^^ . This area denotes the cross section for lensing by a SIS 
lens. Because the one dimensional central velocity dispersion, 



The critical radius her can be expressed equally using either angular 
or proper distances (Kochanek 1993a). 
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o 
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L 
0.0 



0.5 1.0 

separation in arcsecond 



Fig. 3. The ASF for NOT observations at 0'.'7 and I'.'O seeing (FWHM) 
(solid hues). Also included, for comparison, are the ESO-KP ASFs 
for ground-based observations with inspection by EYE (dotted) and 
PSF-subtraction (dashed) for 1.0" seeing and the ASF for Hubble 
Space Telescope (HST) snapshot survey (dash-dot) (see Surdej et al. 
1993). 



(7, is the only free parameter in this model, the statistical analy- 
sis is very sensitive to a. Recently, Kochanek (1993c) proposed 
that the correction factor of ^3/2 to the a value, which is com- 
monly applied in lens statistics to account for the presence of 
DM in late-type galaxies (TOG, F&T), should not be used. We 
will discuss this in Section 6 and present results for both cases 
(with and without the DM correction). 

As usual, we estimate the population of SIS galaxies from 
the Schechter galaxy luminosity function 



L\ L 
jgdL = n^. [ — ] exp(-—)dL/L^, 



(2) 



where the galaxy number density n* = 1.56(±0.4) x 10 ^ 
Mpc-3 and a = -1.1 ± 0.1, following Efstathiou et al. (1988). 
The Hubble constant is Hq= 100 h km s~ ' Mpc~ ' . The galaxy 
luminosity, L, is estimated by using the Faber- Jackson relation 



1/7 



(3) 



where the exponent is taken to be 7 = 4 (Tully & Fisher 1977; 
Faber & Jackson 1976). The value of the exponent 7 is chosen 
(rather than using the tabulated values of 3.7, 3.7 and 2.6 for 
E, SO and S galaxies, respectively) for computational reasons 
and do not affect significantly the statistical results (Kochanek 
1993c). Furthermore, we assume that the number of lenses 
is constant in a co-moving (i.e. non-evolving) volume, hence 
the density n; of lenses is given by n; = no(l + zi)^, where 
the (1 + zi)^ factor arises from the fact that the co-moving 
volume decreases with increasing z, due to the expansion of the 
universe. The present density of isolated galaxies, no, is found 
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Table 3. The relative abundance, central velocity dispersions and computed F values of late-type E/SO galaxies and early-type S spiral galaxies, 
using a. = —1.1 and 7 = 4 (Kochanek 1993c). The values of and uncertainties in the abundance, ra, and cr, are taken from F&T (and references 
therein). Note that the F values are computed using ra, = 1 .56(±0.4) x 10~^ h^ Mpc~' and with/without the DM correction factor of 



type 


abund. 


(7, \kms '] 


F 


(7, \kvas 


-'] 


F 


E 


12% 




0.0085 ± 0.0044 


225%, X 


^3/2 


0.0191 ± 0.0099 


SO 


19% 


206+_'lo 


0.0094 ± 0.0058 


2<d^% X 


072 


0.0212 ±0.0115 


S 


69% 


143+^3 


0.0080 ± 0.0042 


143!.^3 




0.0080 ±0.0042 


E/SO/S 


100% 




0.0259 ± 0.0084 






0.0483 ± 0.0157 


E/SO 


31% 




0.0179 ±0.0073 






0.0403 ± 0.0151 



from integrating the Schecter luminosity function Eq. (2) over 
all luininosities (i.e. no = 4>gdL). 

The optical depth, t, for lensing of a distant QSO with 
magnitude my , and redshift z« , is given by 



(4) 



where F = -^(nocr^) = -^r[l + a+4/'y]n^ai and T is the 
gamma-function (TOG, F&T). 

In Table 3 we list all relevant abundances, velocity dis- 
persions and respective F parameter values with/without DM 
corrections for both E/SO and E/SO/S populations. Note that 
the DM corrections are only applied to late-type galaxies (i.e. 
E/SO). 

In our statistical model we must also take into consideration 
the important magnification bias, which increases the expected 
number of lenses. This is due to the magnification of (strongly) 
lensed images, which effectively alters the magnitude limited 
sample (Narayan & Wallington 1993). We estimate the magni- 
fication bias directly from the QSO luminosity function, which 
we define as in F&T (based on data from Hartwick & Schade 
(1990)) 

r ioO-86(™B-i9.i5) forms < 19.15 
^^(mB)- |jq0.28(™b-19.15) forms > 19.15 " 



(5) 



Here UqimB) represents the counts of QSOs at apparent mag- 
nitude ruB- The differential magnification bias, simply being 
the over-representation of QSOs at magnitude tub due to the 
magnification of originally dimmer QSOs with magnitudes 
ruB + 2.5logM, is given by 



B{mB,M) = 



cLPm n^imB H- 2.5 log M) 



dM 



nJmB) 



(6) 



where dPu/dM = E/M^ is the magnification (M) probability 
distribution for our particular SIS lens model, assuming the 
QSO is a point source. 

5. Results from the NOT sample 

We here present the objectively discovered, possible candidates 
of gravitational lensing from the 168 QSO NOT sample. The 



candidates within the boundaries of the NOT ASF are listed 
in Table 4. Each candidate is also discussed individually and 
additional colour and spectral data, which might enable us to 
reject or approve a candidate, are also presented here. Some 
objects have been observed by other surveys, in which case 
we present their conclusions. Since only a few candidates have 
been observed spectroscopically at NOT using the low dis- 
persion spectrograph (LDS), not many of the candidates have 
been rejected. Furthermore, because no photometric-calibration 
exposures were taken for the NOT sample, we are not able to 
supply colour estimates. We have, however, compared the mag- 
nitude difference between the QSO and nearby companions for 
the various filters and use similar magnitude differences as an 
indication of colour resemblance. 



- PG 1715+535' s companion is confirmed by Yee et al. 

(1993) and Crampton et al. (1992) to be a star, having also 
been investigated spectroscopically by these surveys. Its 
measured colour differences (Amy = 0.12, Amjj = 0.12 
and Amj = —0.05) could not reject either interpretation, 
which is why it could only be determined spectroscopically. 
The small magnitude difference (i.e. high magnification) is 
the main contribution to the high lensing probability. 

- TEX 2127+176 and its companion have a significant colour 
difference (Amy = 1.3, Amjj = 0.5 and Amj = 0.3) and 
is therefore not a strong candidate. Spectroscopic obser- 
vations were unfortunately not properly achieved. Because 
of the small angle separating the QSO (my ~ 18) and 
companion (my ~ 19.5), we suspect that the companion 
was not satisfactorily aligned on the slit. We were able to 
identify the rather strong Civ emission line in the spectra 
of the QSO, while nothing was detectable at that location 
in the fainter spectrum. Further spectroscopic observations 
are planned in 1995 and will hopefully confirm or reject the 
lensing hypothesis. 

- Q 1442+295 and its companion show a significant colour 
difference (Amy = 2.0, Amjj = 1.6 and Amj = 1.1), 
indicating that the companion is quite redder than the QSO. 
Spectroscopic data show no sign of a distinct alikeness in 
emission line or continuum. We conclude that this object 
is not a strong GL candidate, most probably a star. Recent 
spectroscopic observations of this pair, carried out at the 
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Table 4. List of possible GL-candidates in the NOT-sampIe, with red- 
shift, V band magnitude, angular separation and magnitude difference. 
Objects are sorted in decreasing relative calculated probabihty of the 
companion being a lensed image of the QSO. 



Object: 




mv 


e 


A m 


comment 


PG1715+535 


1.93 


16.50 


3'.'7 


0.12 


star 


TEX2127+176 


2.01 


18.00 


2'.'1 


1.30 


red color 


Q1442+295 


2.67 


16.20 


3'.'4 


2.00 


red color 


TEX2 11 6+203 


1.68 


17.00 


3'.'7 


2.75 


red color 


Q1222.7+22.6 


2.29 


18.00 


4'.'3 


0.80 


spectra needed 


4C 56.28 


1.59 


17.30 


3'.'2 


2.70 


red color 


H0905+1507 


3.16 


19.00 


3'.'1 


2.70 


spectra needed 


B3 1621+392 


1.97 


17.50 


3'.'8 


2.30 


star 


Q0308+1902 


2.84 


18.60 


4'.'4 


3.10 


galaxy 


HS0855+2549 


1.80 


18.00 


3'.'3 


4.30 




HS0850+2852 


2.09 


17.50 


3'.'7 


4.50 




HS0727+634 


2.38 


18.50 


4'.'7 


3.80 


spectra needed 



Notes : 

Col. 1: HLQ identification 
Col. 2: Redshift z of the QSO 

Col. 3: V-magnitude of the QSO as tabulated in other QSO 
catalogues 

Col. 4: Angular separation 9 between the QSO and its companion (if 

more than one, chosen as described in Section 2) 

Col. 5: Magnitude difference in the V-band between the QSO and 

companion 

Col. 6: Comment on components nature or characteristics 



NTT, La Silla, Chile by J. Surdej also show little resem- 
blance between the QSO and companion (Surdej 1994). 

- TEX 2116+203 and its companion show a significant colour 
difference (Amy = 2.75, Amjj = 2.2 and Amj = L5). 
Spectroscopic observations give no indication that the com- 
panion is a lensed image of the QSO. 

- Q 1222.7+22.6has only been imaged in the R and V filters. 
The magnitude differences in these colours are, however, 
quite similar (Amy ^ 0.8 and Amjj 0.65) and within 
what we expect to be plausible. The QSO remains, therefore, 
an interesting candidate of lensing. 

- 4C 56.28 and its companion have been observed by other 
surveys, but have not yet been investigated spectroscopi- 
cally. Our photometric data indicate that the companion is 
somewhat redder than the QSO (in agreement with previ- 
ous surveys), but not more than can be expected from a 
reddening effect due to a lensing galaxy (Amy 2.7 and 
AiriR 2.65). We therefore treat this pair as interesting 
and plan to perform spectroscopic observations in 1995. 

- H0905+1507's companion was observed spectroscopically 
by J.Surdej at the NTT and seems to be quite red, though 
no spectroscopic data is yet available (Surdej 1994). 

- B3 1621 +392 and its companion show fairly similar colours 
(Amy K 2.3, AruR ^ 2.25 and Amj ^ 2.4), but the 
companion has been reported by Bahcall et al. (1992) to be 
a G-type star. 



- Q 0308+19 and its companion show fairly similar colours 
(Amy K. 3.1, Amjj k. 3.0 and Amj k. 3.0), but the 
companion has been reported by Crampton et al. (1992) to 
be a galaxy. 

- HS 0855+25 lacks photometric data and demands further 
investigation if its companion's nature is to be revealed. 

- HS 0850+28 and its companion show fairly similar colours 
(Amy f« 4.5, Amjj f« 4.5 and Amj f« 4.5), but has not 
yet been observed spectroscopically. 

- HS 0727+63 and its companion show fairly similar colours 
(Amjj 3.8 and Amj Ri 3.5), but has not yet been ob- 
served spectroscopically. 

6. Cosmological implications 

To make the best possible statistical analysis (i.e. using as large 
a sample as possible) we combine the NOT sample with the 
full optical SCYM-sample presented partially ^ by Surdej et al. 
(1993). Using likelihood functions, we estimate maximum like- 
lihood velocity dispersions, (7* , and lens effectiveness param- 
eters, F, from the combined sample of 784 HLQs for various 
flat cosmologies. In these computations we employ complete 
integral probabilities. This is the probability for a given QSO 
(with apparent magnitude, rrn, and redshift, Zi) to be lensed by 
an intervening isolated galaxy. 

We now wish to express the probability of finding a lensed 
image of a QSO under the constraints imposed by the obser- 
vational conditions. The observational conditions are charac- 
terised by the seeing (FWHM) of point sources in the images. 
The detectable magnitude difference is dependent on the im- 
age angular separation, and is closely related to the seeing. We 
apply the ASF introduced in Section 3, which describes the re- 
lations of detectable magnitude differences, angular separation 
and seeing for our particular survey. Following Surdej et al. 
(1993) and Kochanek (1993c), we can express the probability, 
pf ^, of finding a gravitationally lensed image in the vicinity 
(within rc) of a QSO within the observational window (defined 
by the selection function) by 

SF r \ ScatF „3 , , 

Pi {mi,Zi) = D^^(zi) 

/ Pc{bcr)Bias{mi, M{2bcr), M2)dbcr. (7) 

^0 

Here the critical radius probability is given by pdbcr) = 
= ^ (refer Eq. (4.10) in Kochanek (1993c)) 4. 
Bias{mi, M, M2) is the integral of B{mi,M), given by 
Eq. (6), and the bias integral, using an infinite upper 
limit, is numerically computed by substitution of variables ^. 

^ We have also included, as did Kochanek (1993c), the extended 
HST snapshot sample presented in Maoz et al. (1993). 

We have expressed the lensing probability in terms of the critical 
radius, bcr, instead of the angular separation 9. 
' There was little or no difference in our results using either this 
method or that of setting the upper limit to a large number, i.e. M2 = 
10000 
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Table 5. Estimated values of cr, and F from the frequency of lensing and the configuration likeUhoods of the observations. The results are shown 
in each row for different combinations of galaxy population with and without DM correction. The 90% confidence level results are shown in 
the upper part of the table for the Einstein-de Sitter universe model (Qm = 1). In the lower part, we give the 90% confidence level results 
for the maximum likelihood value of the density parameter (i.e. Qm = 0.95 for E/SO/S and Qm = 0.25 for E/SO populations). The maximum 
Ukehhood value is indicated by J'™"'" and 0-^'"" , respectively. 



Qm galaxy types 


plow pmax pfixgfi 


lotv -max up 


1.00 E/SO/S 

1.00 E/SO/S (without DM) 
1.00 E/SO 

1.00 E/SO (without DM) 


0.02021 0.03735 0.06589 
0.00898 0.01660 0.02928 
0.01490 0.02477 0.04031 
0.00662 0.01101 0.01792 


164.5 191.8 221.1 
134.3 156.6 180.5 

204.3 232.0 262.0 
166.8 189.4 213.9 


0.95 E/SO/S 

0.95 E/SO/S (without DM) 
0.25 E/SO 

0.25 E/SO (without DM) 


0.01964 0.03615 0.06360 
0.00873 0.01607 0.02826 
0.00609 0.01111 0.01944 
0.00271 0.00494 0.00864 


163.3 190.3 219.1 

133.4 155.3 178.9 

163.3 189.9 218.3 

133.4 155.0 178.3 



Note also that Bias(mi, M(2bcr), M2) is normalised so that 
J^^^^j^^dPM I dMdM = 1. The minimum detectable magni- 
fication (i.e. maximum magnitude difference) is given by the 
expression for the magnification M{2bcr) = ^ j(2b°'^)-\ ' where 
filbcr) = f{0) is the flux ratio corresponding to the detectable 
magnitude difference Am(9) (given by Eq. (1)). In our calcu- 
lations, we have taken advantage of the "Numerical Recipes in 
Fortran" implementation of the Romberg method for integra- 
tion (Press et al. 1992). Note also that in absence of any selec- 
tion function the minimum detectable magnification would be 
2.5Iog2 (TOG, Surdej et al. 1993). 

The probability for finding a lensed pair with a certain red- 
shift and magnitude difference within the ASF, is thus the proba- 
bility for lensing weighted by the ASF (Eq. ( 1 )) for each critical 
radius The configuration probability pdbcr) (the relative 
probability of a separation) multiplied by the ratio of the bias 
including the ASF to that with no ASF, is thus integrated over 
the survey area (r < Tc) and multiplied by the general probabil- 
ity of lensing, leading to Eq. (7). The sum of the probabilities of 
finding a GL-system within the detectable ranges of the selec- 
tion function for each observed object, ^ pf^, will then give 

i 

US the expected number of lensed QSOs in the sample. 

The combined NOTh-SCYM sample, consisting of 784 
QSOs, has been applied in the computations below. Among 
these QSOs, four have been identified as being gravitationally 
lensed: the two double imaged lensed systems "Q 1208H-10H" 
(Magain et al. 1988; Maoz et al. 1992) and "UM 673" (Surdej 
et al. 1987), the clover-leaf "H1413H-117" (Magain et al. 1988) 
and the triple image "PG 1115-H080" (Weymann et al. 1980), 
which were all found in the SCYM-sample. We therefore have 
Nl = 4. The well known "QSO 0957H-561" (Walsh et al. 1979) 
is not within our cutoff radius and is therefore not among the 
GL-candidates for the survey. It is also a well known fact that 
the gravitational potential of this GL-system is not solely due 
to the massive central galaxy, but also due to the surrounding 
cluster (contributing to the large separation of 6.1 arc-second). 
The four lensed systems are included, because even though our 



lens model is not able to reproduce the image configuration of 
two of them, the frequency of lensing is still preserved. 

Kochanek (1993c) estimated (7* for E and SO galaxy popu- 
lations without the DM correction factor, while TOG and F&T 
computed the lens effectiveness parameter, F with the DM cor- 
rection factor. Surdej et al. (1993) estimated the F parameter 
without a priori corrections for DM. Since we do not know 
the importance of S spiral galaxies nor that of dark matter, we 
estimate a■^, for bo th po pulations and with/without the DM cor- 
rection factor of \/3/2. Furthermore, we wish to investigate if it 
is possible to obtain reasonable constraints on the cosmological 
density parameter, Am. 

In these computations we apply the configuration probabil- 
ity likelihood function introduced in Kochanek (1993c) 

Nu Nl Nl 

Ptot=i[ii-pr)l[pri[pck (8) 

i=l i=l k=l 

which, converted into logarithmic likelihoods, becomes 

Nu Nl Nl 

InPtot = -J^pf^ ^^P^ + ^^P'^k (9) 
i=i i=i k=i 

under the assumption that p^^ < 1. It is here assumed that 
Pi , pj and pck are the probabilities of a QSO being not lensed, 
lensed or having a separation b if lensed, respectively. The 
configuration probability of a certain separation b is given by 

T^,{a^,)D'l{Zg) 

Pck = Pc(bcr) 5p Bias{mu , M{bcr), M2), (10) 

Pk 

while pf^ and p^^ are given by Eq. (7). The logarithmic likeli- 
hood function adds together all probabilities and counts lensed 
HLQs as positive and unlensed HLQs as negative contributions. 

We, hence, compute likelihoods for various values of the 
velocity dispersion, a^,, and flat cosmological models (i.e. dif- 
ferent Am = 1 — ^^A values). The value of Am affects the 
proper distance, which is part of the optical depth (Eq. (4)), 



A.O. Jaunsen et al: The NOT GL survey for multiply imaged quasars 



11 




Fig. 4. Likelihood contour plots for various flat cosmologies {Um = 1 — S^a) and velocity dispersions (cr,), assuming lensing is significant for 
all lens types (E/SO/S) (left) and only for late type galaxies (E/SO) (right). Both the NOT+SCYM and SCYM samples are represented m the 
plots by soUd and dotted hnes, respectively. The maximum UkeMhood model for these two samples are indicated by * and +, respectively. 



while (7* governs the F parameter (the optical depth t) and the 
image separation distribution, pdbcr)- The individual lensing 
probability pf^ is computed for each object for 1600 various 
combinations of Am and (7* (a 40 x 40 array). 

We present the likelihood results from these computations 
as contours of the 68%, 90%, 95% and 99% confidence levels 
in Fig. 4. These two-dimensional confidence levels are simply 
the corresponding inverse percentage levels of the maximum 
likelihood. The computation of confidence intervals for the (7* 
or F parameter (a single parameter), can be determined from the 
distribution with one degree of freedom. The 90% confidence 
intervals are therefore determined from the 26% peak values of 
the two dimensional likelihood contours in Fig. 4 (see §5 in 
Kochanek (1993c) or Lupton (1993), p.70). 

In Fig. 4 we have plotted the various confidence levels for 
the NOT+SCYM sample as well as those for just the SCYM 
sample. Note that all of the identified lenses were discovered in 
the SCYM sample, so there is no loss in the number of lenses 
when reducing the sample from NOT+SCYM to SCYM. There 
is, however, a slight difference in the results, specifically (and 
not surprisingly) the uncertainties are somewhat lower for the 
NOT+SCYM sample. Also, because there is a smaller number 
of lenses pr. QSO in the NOT+SCYM sample, the contours 
show that the typical a■^, values are somewhat lower than for the 
SCYM sample. 

We present, in Table 5, the confidence intervals for the F 
parameter and a^, in the case of lensing by all galaxy types 
(E/SO/S) and for late-type galaxies (E/SO). Our results are in 
best agreement with the models (in Table 3) excluding the DM 
correction factor. The results are, however, too uncertain to 
draw any definite conclusions about the DM correction factor. 
If, however, the DM correction factor is excluded, the effect of 



S galaxies is not negligible and should be included as part of 
the lensing population (see Section 4). 

If we inspect the diagrams of Fig. 4, we realize that trying 
to constrain the cosmological models (CIm), will still give too 
uncertain conclusions. The 90% confidence levels on a■^, (with- 
out DM) from our results, however, give 164.5 < (7* < 221.1 
kms"' (for Om = 1-0) and 204.3 < (7* < 262.0 kms"' (for 
Am = 1-0) for E/SO/S and E/SO galaxy populations, respec- 
tively. 

To investigate roughly how the uncertainties are decreasing 
with an increased number of QSOs, we have made simulations 
using synthetic samples by simply increasing the present num- 
ber of QSOs and lenses by a factor, n, of two, five and ten, 
albeit this obviously introduces some unwanted selection ef- 
fects. These tests show that the uncertainties on a■^, diminish 
approximately proportionally to ^/n, while the uncertainties in 
CIm diminishes somewhat more rapidly. 

Conclusively we find it difficult to discriminate between 
various flat cosmologies based on the present statistical model 
and sample, while the possible values of galaxy velocity disper- 
sions, (7*, are more restricted. Further examinations of enlarged 
samples of HLQs and QSOs in addition to improved statistical 
analysis will definitely better discriminate between the various 
parameters. 
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